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SECTION  1 
INTRODUCTION 


Analytic  models  which  represent  photoelectron  currents  in  SGEMP 
calculations,  are  useful  for  overcoming  grid  limitations  of  three-dimensional 
simulations  and  also  to  circumvent  the  large  cost  of  particle  calculations 
in  three  or  two  dimensions.  The  P driver,  an  analytic  model  which  represents 
the  buildup  of  the  space-charge-limited  boundary  layer,  has  been  used  quite 
successfully  in  some  code  calculations.1’2  In  other  calculations  it  has  not 
been  completely  successful.3’4  This  note  examines  the  applicability  of  the 
usual  P driver,  for  SGEMP  calculations,  under  two  limiting  circumstances. 

A source  current  is  prescribed  in  the  first  spatial  grid  above  the 
emitting  surface  in  the  P driver  model.  This  source  current  represents  the 
total  integrated  strength  of  the  source  currents  above  the  emitting  surface. 
If  photoelectrons  travel  too  far  from  the  surface  during  times  of  interest  or 
travel  very  little,  relative  to  grid  size  in  the  code,  the  P driver  approxi- 
mations are  violated.  This  note  examines  the  applicability  of  the  usual  P 
driver  code  model  when  Maxwel 1- solver  spatial  grid  dimensions.  Ax,  are  larger 
than  the  characteristic  dimensions  of  the  SGEMP  boundary  layer,  A^.  In  this 
limit  the  source  region  is  artificially  extended  above  the  surface.  An 
examination  of  the  model  is  also  made  under  the  circumstances  that  a signifi- 
cant number  of  particles  escape  to  distances  of  the  order  of  satellite 
dimensions  during  times  of  interest. 

The  examinations  are  made  by  means  of  two  quasi-static  analytic 
solutions  for  fields  on  the  surface  of  a conducting  sphere.  The  sphere  is 
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chosen  because  the  analytic  solutions  are  available  or  easily  obtained;  the 

quasi-static  limit  is  applicable  for  satellites  which  are  highly  damped 

electromagnetical ly.  The  results  of  the  investigation  show  that  if  Ax  is  a 

large  fraction  of  characteristic  body  dimensions,  R (R  is  the  radius  of  the 

sphere,  for  example),  while  is  a small  fraction  of  R,  the  skin  currents  are 

in  error  to  first  order  by  approximately  Ax/R.  These  results  suggest,  for 

circumstances  where  Ax  is  a significant  fraction  of  body  dimensions  (but 
2 

(Ax/R)  can  be  neglected) , that  the  skin  currents  generated  from  a P calcula- 
tion be  augmented  by  a factor  like  1 + BAx,  where  B is,  in  general, 
dependent  upon  position  (B  = R for  a sphere,  far  from  the  source  region 

boundary).  Unfortunately  for  bodies  other  than  a sphere,  6 is  not  well  defined 

2 

and  most  bodies  have  more  than  one  characteristic  dimension.  If  (Ax/R) 
cannot  be  neglected  there  is  probably  no  simple  correction  factor  available. 

If  in  the  real  physical  problem  a large  fraction  of  particles 
actually  escape  the  boundary  layer,  to  distances  comparable  to  the  satellite 
dimensions,  during  times  of  interest,  their  consideration  may  significantly 
affect  the  skin  current.  Modeling  these  distant  currents  with  a P driver 
alone  forces  them  to  contribute  to  the  skin  currents  as  if  they  were  right 
at  the  surface,  rather  than  at  a distance  away  from  the  surface,  it  is  shown,  in 
this  note,  that  the  error,  in  using  a P driver  to  model  distant  photoelectrons,  is 
largest  at  points  on  the  surface  near  the  source  region  boundary;  the  calculated 
currents  will  be  too  large  if  an  accurate  P calculation  is  made  which  includes 
all  the  geometrical  effects.  (These  effects  reduce  the  fields  which  retard 
photoelectron  motion  and  therefore  increase  the  P above  that  calculated  from 
purely  one-dimensional  considerations.)  To  correct  for  the  effect  of  distant 
electrons  prescribed  current  sources  should  be  defined  in  a portion  of  the 
space  outside  the  exposed  surface,  and  not  just  at  the  exposed  surface.  A 
one-dimensional  code  calculation  which  includes  gometrical  effects  may  still 
be  useful  to  define  these  additional  prescribed  sources.  The  one-dimensional 
calculation, alone, will  be  useful  if  edge  effects  can  be  neglected. 
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In  Section  2 we  briefly  describe  the  usual  method  of  incorporating 
a P driver  in  an  SGEMP  code  calculation.  We  then  examine  two  quasi-static 
sphere  calculations  in  Section  3,  for  Ax  > A^.  In  Section  4 we  examine  the 
situation  in  which  the  effect  of  photoelectrons  at  large  distances  must  be 
considered.  The  spatial  extent  over  which  additional  prescribed  sources 
must  be  defined  is  also  considered.  Section  5 is  a summary  of  results. 


SECTION  2 

CODE  MODEL  OF  P DRIVER 


The  theory  of  the  P driver  is  based  upon  the  assumption  that 
during  the  times  of  interest  photoelectrons  emitted  from  a continuous 
electromagnetic  surface  do  not  go  far  from  the  surface.  That  is  if  R is 
the  characteristic  dimensions  of  the  emitting  surface  and  d(t)  is  the 
distance  the  majority  of  electrons  have  gone  in  a time  t,  then  P is  an 
appropriate  driver  if 

R » d (t)  , 

and  we  are  making  magnetic  field  observations  at  points  which  are  distant 
(greater  than  d(t))  from  discontinuities  (corners  of  an  emitting  object, 
for  example)  in  the  emission  surface.  Under  these  circumstances  the  error 
made  in  treating  the  spatial  photoelectron  currents  as  if  they  occur  only 
at  the  emission  surface  is  small  (0(d/R)).  In  a finite  difference  code  a 
current,  J , equivalent  to  the  spatial  distribution  of  currents  is  prescribed 
in  the  first  half  cell  above  the  emission  surface.  The  prescription  is 
made  by  recognizing  that  the  physical  quantity  which  generates  time  varying 
fields,  for  R > d,  is  the  time  rate  of  change  of  dipole  moment  per  unit 
area  P.  The  idea  of  the  code  model  is  to  simulate  the  P of  the  physical 
problem  by  setting 

Ax  Je  = P , (1) 

where  Ax  is  the  grid  size  perpendicular  to  the  emission  surface.  This 
prescription  essentially  establishes  a charge  density,  whose  center  is  one 
Ax  above  the  surface,  equal  to  P(t)/Ax  as  a function  of  time.  P(t)  is 
obtained  by  means  of  a one-dimensional  particle  calculation.  Problems  can 
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arise  with  this  prescription  if  A^R  is  negligible  when  compared  to  one 
but  Ax/R  is  not,  and  also  when  d(t)  ~ R during  times  of  interest.  The 
computational  model  does  not,  then,  correspond  to  the  physical  situation 
In  the  next  section  we  investigate  the  errors  generated  by  this  lack  of 
correspondence. 
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SECTION  3 

ANALYTIC  MODELS 


Taking  a sphere  as  a representative  object  for  SGEMP  analysis, 
we  investigate  the  effect,  on  skin  currents,  of  radially  extended  dipole 
layer  sources  at  the  surface.  We  assume  that  the  actual  dipole  layer  thick- 
ness is  much  smaller  than  the  radius,  R,  of  the  sphere  and  evaluate  the 
error  generated  when  a dipole  layer  is  set  up  in  accordance  with  Equation  1. 
According  to  this  prescription,  a charge  layer  is  set  up  a distance  Ax 
from  the  surface.  It  is  now  assumed  that  Ax/R  is  not  negligible. 

Reference  5 describes  the  exact  quasi-static  analytic  solution  for 
skin  current  generated  on  a perfectly  conducting  sphere  by  spatial  currents. 
This  solution  is  made  in  terms  of  an  expansion  in  Legendre  polynomials.  In 
Section  3.1  we  use  this  solution  to  evaluate  the  error  due  to  a finite  Ax/R. 
It  will  be  shown  that  to  first  order  in  JIAx/R  the  percent  error  in  the  magni- 
tude of  the  Legendre  coefficient  of  skin  current  is  . 5£  pp.  The  error 
increases  with  the  order  of  the  solution  coefficients  and  can  be  large  where 
higher  order  terms  are  important. 

To  obtain  a physical  understanding  of  how  skin  currents  can  be 
smaller  than  expected  we  consider  a second  solution,  the  skin  currents  on 
a sphere  due  to  a single  point  dipole.  This  solution  can  be  obtained  without 
the  use  of  Legendre  polynomials  and  although  the  solution  is  less  general 
it  lends  itself  to  greater  physical  interpretation. 
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Assuming  that  the  spatial  currents,  J(r,0,t),  are  only  in  the  r 
direction,  substitution  of  these  currents  into  Equation  24  of  Reference  5 
gives  the  equation  describing  the  dependence  of  the  skin  current  coefficient 
Kp (t)  on  the  Legendre  coefficients  of  the  spatial  current  -Jp  as 

CO 

K£(t)  = - ~ J J£(r',t)r''\ir'  , A i 1 . (2) 

R 

In  Equation  2,  r is  the  radial  coordinate  of  the  spherical  coordinate  sys- 
tem whose  center  is  at  the  center  of  the  sphere,  (0  is  the  spherical  polar 
coordinate).  The  currents  are  taken  to  be  axisymmetric  so  the  skin  current 
is  expanded  in  Legendre  polynomials  of  order  one  and  .J  is  expanded  in 
Legendre  polynomials  of  order  zero.  The  Legendre  coefficients  for  this 
latter  expansion  are  ,J^.  The  radius  of  the  sphere  is  R. 


We  now  assume  that  J(r,0,t)  can  be  represented  by 
J(r,0,t)  = - P/Axh(0)f(t)  R £ r s Ax 
J(r,0,t)  =0  Ax  < r < » . 


(3) 


Equations  3 are  essentially  saying  that  the  space-charge  layer  sets  up  with 
the  same  time  history  f(t),  for  any  0,  and  retarded  effects  are  unimportant 
(h(0)  is  independent  of  time).  Substituting  Equations  3 into  2,  after 
defining  6 = Ax/R,  we  find  that 

PY.f(t) 

(t:)  = Cp/6)  i a 2 , (4) 

e(M)  = (1  - (1+6)'U1)(6(A-1))'1  , (5a) 


e(A,5)  a 1 - -y-  + A(A+1)  + 0((£6)3)  , 


(5b) 


and  that 


9 


KjCt) 

= Py^ftlE^fi)  , 

(6) 

^(<5) 

= 6-1£n(l+6)  , 

(7a) 

6 62 

= 1 - — + — . 

2 3 

(7b) 

and  6, 

are  the  Legendre  expansion 

coefficients  of  h(0) 

If  6 •>  0 in  Equations  5 and  7 (6 ) 4 1 so  that  P6^f(t)£  is  the  skin  cur- 

rent coefficient  for  a very  thin  dipole  layer  and  e ^ (6 ) is  the  fraction  of 
the  skin  current  coefficient  observed  for  an  extended  dipole.  Table  1 
calculates  (e^(6))  1 for  6 = .3  from  Equations  5a  and  7a.  is  the 

factor  that  the  observed  value  of  skin  current  coefficient  would  be  multiplied 
by  to  obtain  the  current  due  to  a thin  dipole. 


Table  1.  Correction  factor  as  a function  £,  6 = .3 


£ 

1 

2 

3 

4 

5 

6 

1.14 

1.30 

1.47 

1.64 

1.85 

2.04 

Table  1 shows  that  the  error  made  in  extending  a small  dipole  causes  the 
skin  current  to  be  too  small.  It  also  shows  that  at  those  points  where  a 
large  number  of  polynomials  are  necessary  to  describe  the  solution  (edges, 
boundaries,  etc.)  the  error  is  largest. 

3.2  POINT  DIPOLE 

We  now  consider  the  skin  current  on  a perfectly  conducting  sphere 
due  to  negative  point  charge,  -q,  at  0 = 0,  r = R + Ax  H y.  The  physical 
content  of  this  solution  is  clear  and  it  is  general  in  the  sense  that  a sum 
of  point  dipoles  can  yield  a total  solution.  It  is  simplest  to  begin  from 
the  expression  for  the  electric  field,  E(r,t),  in  space,  finding  the  charge 
density  and  then  the  surface  current.  The  electric  field  is 

10 


E = 


q — (r-k  — ) q C 1 -— )r 

q(r-ky)  H y y * Y 


r-ky| ' 


r-k 


R 


2 3 


(8) 


We  are  assuming  in  Equation  8 that  the  charge  -q  was  removed  from  a neutral 
sphere  leaving  a charge  q on  the  sphere.  The  charge  density  0(9, t)  on  the 
surface  is  then  easily  seen  to  be,  if  e = R/y 


o = - q/4iTR"2[(e3-e)  (l+e2-2ecos6)'3/2-l+e] 


(9) 


The  equation  of  continuity  on  the  sphere  is 


1 3 c • OKI  . 9o 

Rsin0  90  (sin6K)  at  ’ 


(10) 


where  K is  the  skin  current.  Substituting  (9)  into  (10),  integrating  and 
requiring  K(0)  = 0 we  find  that  the  current  1(9)  crossing  a circle  made  by 
the  intersection  of  the  plane  z = Rcos0  and  the  sphere  is 

1(0, t)  = 2TrRsin0K(0,t) 

= q/2(l-e)[" ?-1-  + -C j-pf  + cos  el  , (11) 

L(l+e  -2ecos0)  " J 

where  a dot  over  q means  differentiation  with  respect  to  time.  Note  that  at 
9=0,  1(0)  = q or  our  solution  represents  the  situation  in  which  a thin 
current  exists  between  the  charge  atr=R+Ax,  0=0  and  the  surface  of  the 
sphere  at  0 = 0.  The  magnitude  of  this  current  is  q.  Equation  11  will  form 
the  basis  of  the  analysis  of  the  effect  of  distant  electrons,  in  Section  4. 

In  this  section  we  again  look  at  the  situation  in  which  Ax  > 


Since  we  want  our  solution  to  correspond  to  a charge  near  the 
sphere  we  transform  the  independent  variable  of  Equation  11  from  e to  6 where 

6 = Ax/R  , (12) 


and 
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where 


e = (1+6)'1  , 


1(6, t)  = P(2R)~e(6,0) 


e(6,0)  = (l+6)_1[(2+6) (4(l+6)sin20/2+62)-1/2+cos0]  , (15) 


P = qAx 


If  we  assume  that  6 is  small  and  0 > 6,  then 

2 + 6 1 

(4(l+6)sin20/2+62)1/2  sin0/2 

in  Equation  15  and 


(col:  0/2) 


2 2 

e(6 ,0)  * S(0)  (1-6+62)  - t 

sin30/2 


where 


S(0)  = (sin0/2)  A + cos0  . (19) 

It  is  clear  from  Equation  18  that  as  Ax  + 0 (in  such  a way  that  P is  finite) 

the  distribution  of  current  on  the  sphere  is  given  by  I (9 , t ) = (2R)~*PS(0). 

If  6 is  a significant  fraction  then  the  percent  error  due  to  ...  finite  Ax  is, 

2 2 

from  Equation  18,  just  6 itself,  to  0(6  ).  The  6 error  term  .n  Equation  18 
suggests  that  the  largest  percent  error  will  occur  for  0 near  zero  (near 
the  dipole  position).  Table  2 shows  values  of  S(9)/e(0,6),  computed  from 
Equations  15  and  19,  for  6 = .3,  at  various  0,  to  illustrate  the  difference 
between  currents  due  to  a dipole  of  negligible  length  and  those  due  to  a 
dipole  of  finite  length. 


12 


Table  2.  Correction  factor  for  finite  length  dipole,  6 = .3. 


IT 

32 

IT 

16 

IT 

8 

7T 

4 

3tt 

8 

IT 

7 

5tt 

8 

3 

4 17 

h 

IT 

S(8)/e(6,.3) 

3.39 

2.04 

1.51 

1.34 

1.32 

1 .31 

1 .30 

1.31 

1.30 

1 

S(0)/e(0,o)  is  the  factor  the  currents  generated  by  the  finite  dipole  should 
be  multiplied  by  to  obtain  the  current  due  to  dipole  of  negligible  length. 
Table  2 substantiates  the  suggestion  that  the  error  due  to  a finite  dipole 
length  is  largest  near  the  dipole  (8  s Ax/R)  and  about  Ax/R  for  most  of 
the  region  where  Ax/R  >0.  If  the  sphere  were  excited  by  a distribution  of 
dipoles  Ax  could  .correspond  to  a finite  grid  size;  outside  the  source  region 
the  error  due  to  approximating  a negligibly  thin  dipole  layer  by  one  of  finite 
size  would  then  be  Ax/R,  neglecting  terms  of  0(Ax/R)  ). 
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SECTION  4 

PHOTOELECTRONS  AT  LARGE  DISTANCES 


In  this  section  we  investigate  the  errors  involved  in  computing 

skin  currents  by  means  of  a P driver  when  the  condition  R >>  d(t),  discussed 

in  Section  2,  is  not  satisfied.  To  accomplish  the  investigation  we  analyze 

how  currents  in  space,  at  various  distances  from  the  surface  of  a sphere, 

• 3 — 

contribute  to  the  skin  current  on  that  sphere.  Letting  q -+  — a(r,t)dr  in 
at 

Equation  11,  where  a is  a charge  per  unit  length,  we  can,  by  means  cf  the 
equation  of  continuity  and  an  integration  by  parts,  show  that 

00 

1(0, t)  H J J(x,t)  ^~G(x,0)dx  , (20) 

0 


where 


G(x, 0)  = j 6(1+6) 


1 


2 + 6 


(4(l+6)sin 


.2 


S2)1'2 


+ COS0 


)■ 


(21) 


G in  Equation  21  is  just  e(6,0)  in  Equation  15  multiplied  by  1/2;  6 is  now 
x/R,  however,  and  x is  the  distance  of  the  source  current  from  the  surface 
of  the  sphere.  We  can  easily  show,  from  Equation  21,  that 


where 


9G  _ _1_  ^ 
0x  ” 2R 


± r 

2R  L 


(1+6)' 


63(2a-l)  +62(6ct)  +6(12a)  + 8a1 
(l+6)2(62+6(4a)  + 4a)3/2  -> 


(22) 


B = cos©  , (23) 

- -20 

. a = sin  j . (24) 


14 


Figures  1 through  5 are  plots  of  G as  a function  of  6 for  various 
0 from  0 to  t and  indicate  how  effective  spatial  currents  are  in  producing 
surface  currents.  Figures  1 through  5 demonstrate  that  an  infinitesimal 
current  source  close  to  the  surface  of  the  sphere  (6  = 0)  are  much  more  ef- 
fective in  producing  currents  near  the  source  boundary  (the  source  is  at  0 = 0) 
than  far  from  it.  This  is  no  surprise.  When  the  source  current  is  placed 
farther  from  the  sphere  it  becomes  ineffective  in  producing  currents  near  the 
source  but  less  ineffective  in  producing  currents  farther  from  the  source; 
in  Figure  1,  for  example  G at  6 = .5  is  greatest  for  0 = 3tt/15  and  smallest 
for  0 = tt/15.  As  6 increases  the  change  in  effectiveness  of  a source  current 
as  it  is  placed  further  from  the  sphere  decreases.  These  results  are  con- 
sistent with  those  of  Section  3. 


At  any  angle  however  a source  current  is  more  effective  near  the 

surface  than  farther  from  it.  In  other  words  assuming  that  all  the  emission 

• . 00 

current  is  placed  at  6 = 0,  as  in  the  P model  (P  =Q/  Jdx) — whereas  in  fact 
it  is  distributed  in  space — causes  an  overestimate  to  be  made  in  the  magnitude 
of  the  skin  current  at  any  point  on  the  sphere  surface.  Calculations 5 
showing  a P prescribed  current  source  model  giving  too  small  a result  occur 
because  P,  calculated  from  purely  one-dimensional  considerations  is  too  small 
or  edge  effects*are  important.  Reference  3 suggests  that  the  electric  fields, 
in  the  vicinity  of  the  source  region,  which  retard  the  photoelectrons  motion 
away  from  the  surface,  are  reduced  by  skin  currents.  This  reduction  allows 
more  electrons  to  get  farther  from  the  surface  increasing  P.  To  first  order 
the  reduction  in  the  force  on  the  photoelectron  is  given  by  a factor  of  1 - b 
for  this  later  (the  skin  current  effect)  effect  where  b is  a constant.  The 
reduction  of  the  force  on  the  electron  due  to  the  decrease  of  force  with 


D 


* By  edge  effect  we  mean  spatial  currents  which  are  basically  two-dimensional 
in  nature  and  are  caused  by  electrons  escaping  the  system  from  an  edge  or 
by  electrons  which  are  emitted  near  the  edge  of  the  source  region  but 
return  to  an  unilluminated  portion  of  the  satellite.  Edge  effects  are 
not  considered  in  this  note. 
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D 

distance  from  a finite  body  is  also  of  the  order  -jp  Most  likely  neglecting 

both  these  effects  contribute  to  the  error  in  P.  Calculating  P by  means  of  a 

one-dimensional  code  in  spherical  or  cylindrical  coordinates  would  give  a more 

accurate  P.  It  has  also  been  suggested  that  the  skin  current  effect  could 

be  incorporated  into  the  one-dimensional  calculation*, 6 allowing  an  even  more 

accurate  calculation  of  P to  be  available,  in  those  cases  where  electrons 

An 

do  travel  an  appreciable  fraction  of  the  satellite  dimensions  (—  not  small), 

K 

during  the  times  of  interest. 


Figures  1 through  5 indicate,  however,  that  if  particles  do  move  far 
from  the  surface  and  P now  accurately  represents  the  real  P the  skin  currents 
calculated  by  a P model  will  be  too  large.  An  effective  way  to  compensate 
for  the  overestimate  is  to  prescribe  currents  in  the  space  above  the  emission 
surface  as  well  as  at  the  surface.  These  additional  spatial  currents  could 
be  obtained,  just  as  the  original  P model  was  obtained,  by  running  a one 
dimension  calculation.  The  one-dimensional  calculation  would  be  more  sophisti- 
cated only  in  as  much  as  the  geometrical  effects,  as  described  above,  were 
included  in  the  calculation. 


To  obtain  some  idea  of  how  far  in  space  the  currents  should  be 
prescribed.  Equation  20  was  integrated  for  a simple  spatial  variation.  The 
spatial  current,  J,  was  constant  for  0 ^ x s «.  We  denote  the  skin  current 
caused  by  the  spatial  current  from  0 < x < x'  by  I (0 , x • ) , where  from  Equation 
20  and  22  (dropping  the  t) 


1(0, x') 


x ' 

*7 


(2R)  J J (x)  Gdx  , 
0 


(25) 


and  we  also  define 


* B.  Goplen  did  incorporate  the  skin  current  effect  for  the  special  case 
of  a cylinder  (Reference  3),  into  a one-dimensional  calculation.  The 
results  suggest  that  this  effect  may  be  correct  for  the  discrepancy 
between  a particle  pushing  and  a P simulation  of  the  problem  for  the 
radiated  field. 
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(26) 


1(6,6’) 


1(0, x') 

1(0,-) 


where 


6'  = x'/R  . 


(27) 


Figures  6 through  10  are  plots  of  I as  a function  of  for 
various  . It  is  clear  from  the  figures  that  contributions  to  the  current 
close  to  the  source  (small  ) mainly  come  from  currents  close  to  the  surface. 
At  larger  values  of  currents  farther  from  the  surface  make  a stronger 
contribution  to  the  skin  current.  The  physical  situation  resulting  in  Figures 
6 through  10  is  extreme  in  that  the  spatial  currents  extend  all  the  way  to 
infinity  and  only  a very  small  portion  of  the  surface  is  part  of  the  source 
region.  Even  in  this  extreme  case  the  skin  current  would  be  in  error,  by  less 
than  35  percent,  if  we  integrated  only  up  to  one  radii  from  the  sphere,  for 
most  positions  on  the  spherical  surface. 


Some  judgement  is  necessary  in  deciding  to  what  extent  prescribed 
currents  need  to  be  defined  in  space  in  order  to  perform  an  accurate  SGHMP 
simulation.  Equation  25  and  Figures  6 through  10  should  be  helpful  in  mak- 
ing these  decisions.  In  most  circumstances,  where  peak  skin  current  is  of 
primary  interest,  it  will  probably  not  be  necessary  to  prescribe  the  source 
beyond  one  characteristic  radii  of  the  satellite. 
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Figure  10. 


SECTION  5 

SUMMARY  OF  RESULTS 


An  analysis  of  the  limitations  of  the  P driver  was  made  using  a 

quasi-static  analysis  of  a sphere.  Errors  are  caused  by  either  artificially 

over  extending  the  spatial  region  of  the  photoelectron  source  current  or  not 

really  taking  into  account  the  large  extent  of  the  source  current  region. 

Figures  1 through  5 indicate  that  if  a source  is  over  extended  the  calculated 

skin  currents  will  be  too  small,  as  in  the  case  Ax  > XQ  (Ax  is  the  spatial 

grid  size  in  a multidimensional  SGEMP  Maxwell  solver).  In  contrast,  if  a 

source  current  consists  of  distant  electrons  and  they  are  treated  as  though 

they  are  all  near  the  surface  the  calculated  skin  currents  will  be  too 

large.  To  accurately  account  for  the  error  made  when  Ax  > would  require, 

in  general,  an  understanding  of  the  electromagnetic  nature  of  the  satellite 

2 

but  the  error  will  be  of  the  order  Ax/R,  if  (Ax/R)  is  negligible,  where  R 
is  a characteristic  dimension  of  the  satellite.  For  a sphere  the  error  is 
exactly  Ax/R  far  from  the  source  region  boundary. 

When  the  spatial  extent  of  the  affective  source  region  is  large, 
prescribed  sources  must  be  defined  in  the  space  outside  the  satellite  as  well 
as  at  the  surface.  The  spatial  extent,  as  well  as  the  magnitudes  of  the 
sources  can  probably  be  defined  by  running  a one-dimensional  code  (assuming 
edge  effects  are  negligible)  which  incorporates  the  geometrical  effects 
caused  by  the  finite  size  of  the  satellite.  The  sphere  analysis  suggests 
that  a spatial  extent  of  one  satellite  radii  will  probably  be  sufficient 
for  most  situations  in  which  peak  skin  currents  are  sought. 
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